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ABSTRACT 

We have investigated the relation between the orbital periods (Porb) and the 
spin periods (Pg) of wind-fed X-ray pulsars in high-mass X-ray binaries (HMXBs), 
based on population synthesis calculations of the spin evolution of neutron stars 
during the pre-HMXB stage. We show that most of the neutron stars either have 
steady accretion or still reside in the radio pulsar phase when the donor star 
starts evolving off the main sequence. In the former case the values of Ps can 
be decelerated to be ~ 10^ — 10^ s depending on Porb- We briefly discuss the 
possible origin of the Pg — Porb correlation in Be/X-ray binaries, and the existence 
of HMXBs with main sequence donors. We have also investigated the evolution 
of the inclination angle between the magnetic and spin axes of neutron stars in 
a massive binary, suggesting secular alignment of the magnetic and spin axes 
during their evolution. 

Subject headings: stars: early-type — X-rays: binaries — stars: evolution — 
stars: neutron 



1. INTRODUCTION 

High-mass X-ray binaries (HMXBs) can be roughly divided into two types: a supergiant 
star and a Be star are contained respectively. The X-ray source often is a pulsar and powered 
by accretion of the material offered by the companion star. In the supergiant systems, either 
Roche-lobe or stellar wind accretion occurs, while in the Be systems commonly only the 
latter process takes place since the Be star is well inside its Roche lobe (Tauris & van den 
Heuvel 2005). 

The relation between the spin periods Pg and the orbital periods Porb of HMXBs was 
first studied by Corbet (1984, 1985, 1986), who pointed out that there may be a Pg oc P^^/J 
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correlation for neutron stars in systems with radially expanding winds, while Pg oc P^^^ for 
Be systems, van den Heuvel & Rappaport (1987) noticed that the spin and orbital periods 
agree reasonably well with the former correlation for the supergiant systems in which the 
X-ray source is powered by stellar wind accretion. The Pg — -Porb relation has been studied by 
Stella et al. (1986) in terms of the equilibrium spin period, where the corotation radius equals 
to the magnetospheric radius defined by equating the ram pressure of accreting matter with 
the magnetic pressure of the dipole field of the neutron star (e.g., Davidson & Ostriker 1973; 
Lamb, Pethick & Pines 1973). However, in order to obtain quantitative agreement between 
the predicted and observed spin periods for given orbital periods, the assumption of the 
equilibrium spin periods requires a mass accretion rate two orders of magnitude lower than 
observed for supergiant systems (Stella ct al. 1986). Waters & van Kcrkwijk (1989) presented 
a comprehensive explanation for the Pg — Porb relation for HMXBs. For supergiant systems 
they suggested that the difference between the required and the observed mass accretion 
rates can be explained by assuming that the present-day periods are the equilibrium spin 
periods for the stellar wind of the optical star when it was still on main sequence (MS). The 
Pg — Porb correlation for Be/X-ray binaries was suggested to result from the equilibrium spin 
realized by interaction between the neutron star magnetic field and the Be star's equatorial 
wind - the wider orbits expose to the neutron stars to a lower wind density on average and 
hence lower accretion rate. The latter point was criticized by King (1991), who proposed that 
the observed Pg — Porb relation was probably the result of an earlier evolutionary stage, as 
for supergiant binaries. By calculating the angular momentum carried by the Be star's wind 
material, Li & van den Heuvel (1996) showed that for Be/X-ray binaries in relatively narrow 
orbits (Porb ^ 100 days), the equilibrium spin could be reached by angular momentum 
transfer via a disk, formed in the equatorial wind of the Be stars; in wider binaries, the low 
angular momentum of the wind matter prevents the existence of an accretion disk, and the 
neutron star's spin remains close to its previous equilibrium value at an earlier evolutionary 
stage. Zhang, Li & Wang (2004) have calculated the detailed spin evolution of a wind-fed 
neutron star in OB binaries prior to the HMXB phase, by simulating the time-development 
of the mass loss rate and radius expansion of a 20 Mq donor star. 

The above-mentioned investigations are usually either semi-analytic or numerical only 
for individual cases. To better understand the distribution of HMXBs in the Pg — Porb dia- 
gram, an evolutionary population synthesis incorporated with neutron star spin evolution is 
needed. In the present paper, we have conducted a Monte Carlo population synthesis study 
of the spin evolution of a neutron star in a massive binary. Due to the large theoretical 
uncertainties on the formation and evolution of HMXBs, our work is limited to the spin 
evolution of neutron stars before the HMXB phase. The theoretical considerations are de- 
scribed in §2. The calculated results are presented in §3 and their implications for accretion 
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processes in HMXBs are discussed in §4. 



2. MODEL 

2.1. Spin evolution 

We consider a IAMq magnetized neutron star in binaries with a massive MS companion 
star. Wc adopt a simplified version of the theoretical model outlined by Davies & Pringle 
(1981) to calculate the neutron star's spin evolution before steady accretion occurs, as briefly 
described below. 

Case a: the pulsar phase 

The newborn neutron star first appears as a rapidly rotating, radio pulsar after its 
birth in a supernova explosion, provided that its radiation is strong enough to expel the 
wind material coming from the companion star outside the radius of the light cylinder ric = 

cPs/2t[ or the Bondi accretion radius tq = 2GM/v'^ (Bondi & Hoyle 1944), where G is the 
gravitational constant, M the mass of the neutron star and v^o = lO^Vscras^^ the relative 
wind velocity at the neutron star's orbit. The neutron star spins down due to magnetic 
dipole radiation and/or particle emission, 

where / is the moment of inertia, /i — IO^^/isq Gcm^ the magnetic dipole moment, and Qg 
the angular velocity of the neutron star, respectively. 

The pulsar phase will break down in two distinct circumstances. 

(i) The pulsar phase will cease when the wind material penetrates inside ric, and the 
resulting is derived by balancing radiation pressure from the pulsar with the stellar wind 
ram pressure at ric (Davies & Pringle 1981), 

P. 0.84^M,-^/^(M/Mo) V3^-V6 3^ (2) 

where M — lO^^Mis gs~^ is the mass fiow rate onto the neutron star. 

(ii) The pulsar phase will end when the pressure gradient become important as the outer 
boundary i?a of the envelope where the density p ~ poo approaching to tq, i.e., Ra — tq, 
and the corresponding Pg is (Davies & Pringle 1981), 

PaC^1.2M[,'/'^ll^'v-'/'s. (3) 
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Case b: the propeller phase 

The propeller phase follows the pulsar phase. During this period the magnetospheric 
radius = [fi^/{2GMM^)]^/\ is larger than the corotation radius = {GM/fliy/^. The 
plasma interacts with the neutron star magnetosphere, but further accretion is inhibited 
by the centrifugal barrier, and the infalling matter is accelerated outward, taking away the 
angular momentum of the neutron star. Though the propeller effect has be investigated 
extensively, there exist large uncertainties in the efficiency of angular momentum loss during 
the propeller regime (e.g. Pringle & Rees 1972; lllarionov & Sunyaev 1975; Davies & Pringle 
1981; Wang & Robertson 1985; Ikhsanov 2001). Here we assume that the infalling material is 
ejected with the corotating velocity at i?m, and the spin-down torque is (Wang & Robertson 
1985; Jiang & Li 2005) 

N^m,^ -MRlQ,. (4) 
The typical spin-down time-scale r = |r2s/fis| can be estimated to be 

T ~ 2.2 X lOVao'^'^if W^0)'^'^45 yr. (5) 

The spin-down process ends until Pg reaches the equilibrium spin period 

Pec ^ 17/.^ M-'/^(M/Mo)-^/^ s, (6) 

and we assume that steady accretion takes place afterwards (see however, Arons & Lea 1976; 
Eisner & Lamb 1976). 

The spin period of the neutron star may be further changed after P > Pgq in the 
following accretion phase. But we stop the calculations when either Ppq is reached within the 
MS lifetime or the companion star evolves off the MS (so that strong wind accretion or Roche 
lobe overflow occurs) (see also Waters & van Kerkwijk 1989). In this work we do not consider 
narrow HMXBs with Roche-lobe overflow (hke SMC X-1, Cen X-3, and LMC X-4), where 
the neutron stars are most likely to accrete from an accretion disk. For wind-fed systems like 
Vela X— 1, numerical calculations (e.g. Pryxell & Taam 1988; Matsuda et al. 1992; Anzer & 
Borner 1995; Ruffert 1999) have suggested that there are no significant angular momentum 
transfer onto the neutron star when radially-expanding wind material accretes onto the 
neutron star. This may result in only small deviation from the instantaneous (equilibrium if 
reached) spin periods at the beginning of the accretion phase. CGRO/BATSE observations 
have shown a random walk in their spin frequencies with alternating spin-up and spin-down 
(Bildsten et al. 1997). The spin evolution in Be/X-ray pulsars has a completely different 
picture, and will be discussed separately. 
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2.2. Evolution of the inclination angle 



Besides the spin evolution, we have considered the evolution of the inclination angle x 
between the spin and magnetic axes in wind-fed neutron stars. This part of work is partially 
motivated by Bulik et al. (2003). By analyzing the light curves of 89 accretion-powered 
pulsars these authors obtained an upper limit on the inclination angles x < 50°, suggesting 
that the magnetic axis tended to be aligned with the rotation axis. Similar analyzes and 
conclusions have also been conducted by Wang & Welter (1981), Leahy (1991), and Leahy & 
Li (1995). More recently Blay et al. (2005) presented the contemporaneous high-energy and 
radio observations of the HMXB 4U 2206-^54 conducted with INTEGRAL and VLA, which 
firmly indicate that 4U 2206+54 hosts a magnetic accreting neutron star. The absence of 
pulsations are most likely to be due to geometrical effect that the inclination angle of the 
neutron star is very small so that no pulsation can be detected. 

Wang & Robnik (1982) show that the inclination angle in a binary X-ray pulsar will 
vary as a result of its interaction with the surrounding plasma, becoming smaller (larger) 
during spin-down (spin- up). According to Eqs. (27), (36c), and (37c) in Wang & Robnik 
(1982), we can derive the changing rate of x follows. 



where 9c is the critical value of the meridional angle measured from the north magnetic pole. 



We employed an evolutionary population synthesis method to investigate the distribu- 
tion of the orbital periods and companion masses for the natal neutron star binaries (i.e. at 
the moment when the neutron star was born) . We started with a large set of primordial bina- 
ries and generated the systems that evolved to contain a neutron star and a MS companion. 
The initial mass function of Kroupa, Tout & Gilmore (1993) was taken for the masses of the 
primary star (the progenitor of the neutron star of mass Mi). For the secondary star (of mass 
M2), we assumed a uniform distribution of the mass ratio for < g = M2/M1 < 1. A uni- 
form distribution of In a was also taken for the binary separation a. We adopted the binary 
population synthesis code developed by Hurley et al. (2000, 2002) to evolve the primordial 
binaries. This code incorporates evolution of single stars with binary-star interactions, such 
as mass transfer, mass accretion, common-envelope (CE) evolution, collisions, supernova 
kicks, tidal friction and angular momentum loss mechanisms. Most of our adopted parame- 
ters are the same as the standard ones described in Hurley et al. (2002). The star formation 




37Vsinx/(/asin2^e), X < Oc 
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Evolution of the mass flow rate onto neutron star 
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rate parameter is 5" = 7.6085 yr"^, corresponding to a rate of ~ 0.02 yr~^ for core-collapse 
supernovae (SNe) in our Galaxy, assume all the stars with masses > SM© die through SNe. 
During the SN explosions, we apply a Maxwellian distribution in the kick velocities with a 
mean of 265 kms^^ imparted on the newborn neutron stars (Hobbs et al. 2005). The treat- 
ment of Roche-lobe overflow (RLOF) mass transfer in the primordial binary is presented 
specifically in Hurley et al. (2002), and here we describe the stability criterion of mass trans- 
fer briefly. Mass transfer via Roche-lobe overflow takes place on either nuclear, thermal, or 
dynamical time-scales, depending on whether the primary remains in thermal equilibrium as 
it loses mass, and the radius of the primary increases faster than the Roche-lobe. Stars with 
deep surface convective zones, e.g. giants or naked helium giants, are generally unstable to 
dynamical timcscale mass loss and will enter a CE evolution. The stable mass accretion rate 
of the secondary star is limited by Eddington accretion rate. Actually, the secondary may be 
spun up and become a Be star when it accretes enough mass. This circumstance however, 
is not included in our consideration since the origin of Be phenomena is still unclear and it 
is difficult to model the mass transfer processes in Be/X-ray binaries. The CE parameter a 
was set to be 1 as a typical value and we varied it from 0.1 to 2 in the calculations (Dewi & 
Tauris 2000; Tauris & Dewi 2001). Our product is a set of "incipient" neutron star in mas- 
sive binaries with a particular distribution of Porb and Mi (shown in Fig. 1). Other binary 
parameters. sTich as the radii, surface temperatures and luminosities of the companion stars 
can also be obtained. These parameters were then used to evaluate the mass loss rates from 
the companion stars and the mass flow rates onto the neutron stars. 

For comparison with the observed properties of HMXBs, our calculations are hmited 
to systems with Porb < 1000 days and \^Mq ^ M2 ^ 30Mq. The mass loss rate M2 was 
estimated with the prescription described by Nieuwenhuijzen & de Jager (1990), 

-M2 = 9.6 X 10-'^P^-^'L^-24M°-'6 MqYt-\ (8) 

where R2 and L2 arc the radius and luminosity of the companion star. All the quantities in 
Eq. (8) are evaluated in solar units. Assume that the stellar wind expands isotropically at a 
speed of v^, the wind density at the orbit of the neutron star is 

= -M2/(47raV), (9) 
and the mass flow rate onto neutron star is roughly given by (Bondi & Hoyle 1944) 



(10) 
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3. RESULTS 

We have done calculations of the spin and inclination evolutions for 6 x 10^ neutron 
star binary systems based on the theoretical model presented in §2. For the initial neutron 
star magnetic fields B we assumed that log-B is distributed normally with a mean of 12.5 
and a standard deviation of 0.3. No field decay is considered. The initial distribution of the 
inclination angle x is randomly distributed in the range of (0, 7r/2). We set the initial spin 
periods to be distributed uniformly between 10 and 100 ms, though they have httle influence 
on the final results. We stop our calculations when either Pg reaches Peq or the companion 
star begins to evolve off the MS with the arguments afore. 

The calculated results are presented in Figs. 2—5. Panels (a) and (b) in Figs. 2—3 
correspond to the relative wind velocity at the neutron star's orbit vs — 1 and 2 respectively. 
Figure 2 shows the final distribution of all the binaries in the Pg — -Porb diagram. The relative 
numbers of binary systems are displayed sketchily with the darkness of the shaded regions. 
According to our calculation, we find that, when f § = 1, around 68% neutron stars can reach 
the equilibrium periods to allow wind accretion within the MS lifetime ^ms i while others 
would still be in either the pulsar (~ 31%) or the propeller (~ 1%) phases. This can be 
clearly seen in Fig 2, where the former two groups of neutron stars occupy the distinct upper 
and lower shaded regions. The small number of stars in the propeller phase originates from 
its much shorter duration (Davies & Pringle 1981). If Vs is increased to be 2, the mass flow 
rates onto the neutron stars will be lower by a factor of ~ 16 according to Eqs. (9) and (10), 
further extending the spin-down time in the pulsar phase. So the corresponding numbers 
become ~ 42%, 55% and 3%, respectively. The larger wind velocity also induces longer 
equilibrium periods. Figure 2 reveals that up to half of neutron stars in binary systems 
cannot reach the equilibrium periods, indicating that there might be hundreds of "sleeping" 
neutron star + MS star binaries hke PSR B1259+63 in the Galaxy (see also Fig. 4 below). 
They could be observed as radio pulsars or in X-rays depending on whether the interaction 
between the neutron star and the companion's winds becoming active. Most of them have 
Pg ~ 1 s and Porb from tens of to hundreds of days. 

To compare the calculated results with observations of HMXBs, we show the distribution 
of those neutron star binaries with Pg > Peq and of the observed HMXBs in Fig. 3. Asterisks 
and diamonds mark the Be and supergiant wind-fed HMXBs respectively, and crosses for 
Roche- lobe overflow systems (data are taken from Raguzova & Popov 2005). We flrst discuss 
the persistent, supergiant systems. Obviously the spin periods of neutron stars can approach 
Peq if the total spin-down timescale tspin is less than the MS lifetime tus of the companion 
star. As we argued before, for this type of systems, the current spin periods are likely 
to be around the equilibrium periods Peq (MS) attained during the MS stage because of 
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the inefficient angular momentum transfer in tfie accretion pliase. For tfiese systems to 
be observed as HMXBs, tlie accretion rate must liave been increased to be liigli enougli 
tliat tlie spin periods are larger than the (current) equilibrium period Pgq (sg) when the 
companion star becomes a supergiant (Stella et al. 1986; Waters & van Kcrkwijk 1989), which 
is plotted with the dotted line in Fig. 3 according to Eq. (6). Here we take M2 = 20Mq, 
M2 = 10^^ Mqjt^^, and /xso = 3 (Coburn et al. 2004) as typical values of the supergiant 
systems. This period roughly serves as the lower limit for Pg of the neutron stars in supergiant 
HMXBs. Moreover, since the mass accretion rates generally decrease with Porb, only those 
in narrow orbits (Porb i$ 15 d) can have X-ray luminosities > 10^^ ergs~^ as observed. 

Be/X-ray binaries generally have much longer orbital periods than the supergiant sys- 
tems, and are usually transient X-ray sources observable during outbursts. The structure of 
the Be star's winds is quite complicated, consisting of a relatively dense, slowly expanding, 
disk-like, equatorial wind and a fast, isotropically expanding, polar wind (Waters & van 
Kerkwijk 1989). Waters & van Kerkwijk (1989) suggested that the P, — Porb correlation 
for Be/X-ray binaries could be explained by setting Pg = Pcq obtained in the disk winds. 
However, current understanding of the Be disk winds (Okazaki & Negueruela 2001) is very 
different from the radial wind model used by Waters & van Kerkwijk (1989). The evolu- 
tion of the Be star's disk and the high eccentricities of the orbits also mean that the mass 
flow rate onto neutron stars in Be systems are always highly variable on both long-term 
and short-term timescales. Because of these complications we did not calculated the spin 
evolution of neutron stars accreting from the disk winds, but present constraints on their 
possible locations in the Pg — Porb diagram instead. We first note that the calculated Peq 
(MS), under the assumption of isotropic winds, can be regarded as the upper limits for Pg 
of Be/X-ray binaries, since M is much lower from the polar winds than from the disk winds 
(these values should be taken as the averaged ones since we have ignored the eccentricity of 
the Be star binaries). Then, similar as in the supergiant systems, we can set the lower limit 
of Pg for Be/X-ray binaries to be the equilibrium period Peq (disk) in outbursts when the the 
disk winds dominate accretion. From the observed data of 36 Be/X-ray binaries compiled by 
Raguzova & Popov (2005), we have derived a correlation between the maximum luminosities 
and the orbital periods, 

log(Lx,max/10=^'ergs-i) = 4.53(±0.66) - 1.50(±0.33) log(Porb/d). (11) 

Insert Eq. (11) and /^ao = 3 into Eq. (6) we obtain 

log Peq(disk) ~ -0.29 + 0.64 log(Porb/d), (12) 

which is plotted in the sohd line in Fig. 3. We see that a large fraction of Be/X-ray binaries 
he between Peq (MS) and Peq (d). The pecuhar location of A 0535—669 may be due to its 
relatively low magnetic field. 
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The cumulative luminosity distributions of HMXBs based on our population synthesis 
calculations are shown in Fig. 4 with different values of the parameters a, Vg, and a. The 
results show that change of these parameters do not significantly influence the final outcome, 
which arc roughly compatible with the observed luminosity distributions of HMXBs in our 
Galaxy (Grimm et al. 2002). 

In Fig. 5 we plot the inclination angle distribution for X-ray pulsars in HMXBs. It is 
easily seen that most of the inclination angles are within 1 rad, in general agreement with 
the analysis of the observation data (Bulik et al. 2003). The evolution of the inclination 
angle seems to be insensitive to the values of the relative wind velocity. 

4. DISCUSSION 

We have calculated the spin evolution of neutron stars in a massive binary. The main 
ideas in this work are of two folds. First, to appear as X-ray binaries, the spin periods of 
neutron stars should be longer than the current (instantaneous) equilibrium periods. Second, 
to satisfy this condition, the neutron stars should be spun-down enough during the previous 
propeller phase when the companion stars had relatively weak winds. These two requirements 
set constraints on the spin period distribution in wind-fed X-ray pulsars. Since M depends 
on the orbital separations, this may result in the particular distributions of the supergiant 
and Be HMXBs in the Pg — Porb diagram. The numerical results presented in Figs. 2 and 
3, roughly accounting for the main features of the Pg — Porb distribution of HMXBs, lend 
support to the original suggestions by Corbet (1984), Stella et al. (1986), and Waters & 
van Kerkwijk (1989). Our calculations of the alignment during spin-down imply that most 
HMXBs may have relatively small inclination angles. 

Wc note that our preliminary results are subject of lots of uncertainties and simplified 
treatment adopted. In case of isotropic wind accretion, a considerable fraction of neutron 
stars can reach P^q when the companion star is still on the MS. The detail number is de- 
termined not only by the binary evolutionary processes, but also critically by the propeller 
mechanism. The spin-down torque (Eq. [4]) adopted in this work is among the most efficient 
ones (Jiang & Li 2005, and references therein), and the results should be taken as to be the 
most optimistic. For less efficient spin-down torques ones, the number of binaries in the pro- 
peller phase will obviously increase due to the longer spin-down timescales. The magnitude 
of Peq depends on the mass loss rate, the wind velocity, and the magnetic field strength of 
the neutron star. With the mass loss rates given by Eq. (9) and typical magnetic fields of 
~ 3 X 10^^ G, the neutron star's spin can be decelerated to be ~ 10^ — 10'^ s within tus 
of the secondary star. The very long period Pg = 10^ s of 2S 0114-1-650 (Hall et al. 2000), 
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however, may be explained by a ultra-strong initial magnetic field {B ^ 10^^ G), in order to 
be spun-down efficiently by the propeller effect (Li & van den Heuvel 1999). 

The structure of the Be star's winds is much more complex and variable than in su- 
pergiant systems, and the simple propeller spin-down may not be applicable to this type of 
systems. It is likely that the Pg — -forb distribution and correlation of Be/X-ray binaries result 
from the balance between the spin-up during outbursts and spin-down during quiescence. 

Assume that the spin-up and spin-down torques are MoiGMRa^oV^'^ and —M^R'^^^Qs respec- 
tively, we can derive the following expression for the equilibrium period (see also Menou et 
al. 1999), 

^e. - 27r(§^)-/^(5§a)(^) oc M-'^^^nh (13) 

Here the subscripts o and q denote quantities evaluated during outbursts and quiescence, 
respectively. All the three terms on the right hand side of Eq. (13) are hkely to increase with 
Porb; which might account for the observed Pg — Porb correlation of Be/X-ray binaries. 

Our calculated results also indicate the existence of HMXBs with MS donors in the 
Galaxy. These systems, as the progenitors of supergiant HMXBs, are less luminous than the 
latter (but with similar spin periods) . Figure 4 suggests that there could be a few hundred 
of these sources in the Galaxy with X-ray luminosities ranging from ~ 10^^ to 10^^ ergs~^ (If 
the propeller effect is not considered, the total number of X-ray binaries, most of which are 
Bc/X-ray binaries, can reach a few thousand). The HMXB 4U2206-I-54 is likely to be the 
prototype of this kind of sources (Ribo et al. 2005). A number of papers have also suggestted 
that some neutron stars receive low kick speeds of < 50 kms~^ at birth (Pfahl et al. 2002; 
Podsiadlowski et al. 2004; Dewi, Podsiadlowski & Pols 2005). If all the neutron stars are 
born with such small kicks, our calculations show that there would be about 4 — 5 times 
more X-ray sources produced. 
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Fig. 1. — The initial orbital period distribution of neutron star + massive MS star binaries 
(left) and the incipient mass distribution of the MS companion stars (right). 
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Fig. 2. — The distribution of neutron star binaries in the Pg — Porb diagram when Pg 
or at the age of ^ms- The color bars denote the relative numbers in different regions 
left and right panels correspond to = 1 and 2, respectively. 
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Fig. 3. — The Pg— -Porb distribution of wind-fed HMXBs. The dashed and sohd hnes represent 
the lower hmits of the spin periods for supergiant and Be systems respectively. Asterisks and 
diamonds mark the Be and supergiant wind-fed HMXBs, and crosses for Roche-lobe overflow 
systems, respectively. The left and right panels correspond to vq = 1 and 2, respectively. 
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Fig. 4. — The luminosity distribution of wind-fed HMXBs based on population synthesis 
calculations. Top: We take = 1, a = 190 kms~^ {a = 0.1 and 1 in the left and right 
panels respectively). Bottom: We take a = 2, vs = 2 {a = 190 kms~^ and 265 kms^^ in 
the left and right panels respectively). 



-18- 




Fig. 5. — The relative distribution of the inchnation angles of HMXBs. The left and right 
panels correspond to = 1 and 2, respectively. 



